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ABSTRACT 



Aims. We present an analysis of a peculiar supergiant B-type star (VFTS698/Melnick 2/Parker 1797) in the 30 Doradus 
region of the Large Magellanic Cloud which exhibits characteristics similar to the broad class of B[e] stars. 
Methods. We analyse optical spectra from the VLT-FLAMES survey, together with archival optical and infrared pho- 
tometry and X-ray imaging to characterise the system. 

Results. We find radial velocity variations of around 400kms~^ in the high excitation Siiv, Nlll and Hen spectra, and 
photometric variability of ~0.6mag with a period of 12.7 d. In addition, we detect long-term photometric variations of 
~0.25 mag, which may be due to a longer-term variability with a period of ~400 d. 

Conclusions. We conclude that VFTS698 is likely an interacting binary comprising an early B-type star secondary orbit- 
ing a veiled, more massive companion. Spectral evidence suggests a mid-to-late B-type primary, but this may originate 
from an optically-thick accretion disc directly surrounding the primary. 



Key words, stars: early-type - 
Be - stars: variable: general 



binaries: spectroscopic - stars: peculiar (except chemically peculiar) - stars: emission-line. 



1. Introduction 



Massive stars, whilst int rinsically rare due to the slope of 
the initial mass function (jSalpeteilllQSSf) , have an important 
role in influencing their surroundings. As the progenitors 
of core-collapse supernovae (e.g. ^Smartt 2009), they inject 
both nuclear processed material and energy into the inter- 
stellar medium. Additionally their high luminosity permit 
us to observe them beyoiid the Milky Way ( e.g. iMassevI 
120031: iLevesaue et al.l [20051: iMassev et al.ll2009D . providing 
a crucial test of stellar evolution at different metallicities. 

Variability is a defining characteristic for several of the 
subclasses of massive stars, e.g., luminous blue variables 
(LBVs) and Cepheids. Such variability can take a wide 
variety of forms: pulsations in Cepheids can give rise to 
stable, periodic light-curves, while LBVs undergo irregular 
and unpredi ctable eruptive outbursts (e. g. the most famous. 
Eta Carina, Ivan Genderen fc"Th| [T98?) which can lead to 
sudden increases in magnitude. Individual types of massive 
stars can also display more than one type of variability, for 
example some LBVs can also display gradual changes in 
magnitude over timescales of years. 



Spectroscopic and photometric variability in massive 
stars can also be caused by binarity, with as many as 40% 
of massive stars found to have binary companions (B-type: 
Dunstall et al. to be submi tted, 0-type: Sana et al. to be 
submitted). iTarasovl (200(f) presented a review of short- 
period interacting binaries which give rise to emission-line 
phenomena during the course of their evolution. Short- 
period variability (~5 d) in emission features is associated 
with "classical Algol" objects, whereas variability of the 
order of tens of days can be seen in W Serpentis bina- 
ries or double periodic v ariables (DPVs. IPlavec et aI1ll980l : 
iMennickent et al.ll2003l) . 

B[e] supergiants (sgB[e]) are the most homogeneous and 
distinct group of s tars that show the '[e]' phenomenon 
(jLamers et al.l[T998h . Their spectra are characterised by the 
presence of narrow emission from both permitted and for- 
bidden lines of iron-group elements. SgB[e] stars have high 
luminosities (log L/Lq > 4) and are therefore easily observ- 
able in the Magellanic Clouds. Ten sgB[e]s in the Large 
Magellanic Cloud (LMC) are well documentecQ, allowing 
concise crit eria to be pla.ced ii pon the characterisation of 
these stars (|Lamers et al.lll99"8h . The y show a possi b le evo - 
lutionary connection to LBVs, e.g., iHerrero et all (|2010[ ). 



* Based on observations at the European Southern 
Observatory Very Large Telescope in programme 182.D-0222 



^ iLamers et all (|1998|') listed 11 B[e ] supe rgiants in the LMC, 
later revised to 10 bv lBonanos et all (120091 1. 
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defined V39 in tlie Small Magellanic Cloud (SMC) as ei- 
ther an LBV candidate or a sgB[el. In addition, sgB[cls may 
be the result of binary mergers (jPodsiadlowski ct al. 2006), 
although their exact nature remains unclear at present. 

In this paper we document the spectral and pho- 
tospheric characteristics of an intriguing radial velocity 
(RV) variable from the VLT-FLAMES Tarantula Survey 
(|Evans et all 1201 iL VFTS; hereafter Paper I). Th e target 
VFT S698, has previous ly been know n as Melnick 2 (jMelnickl 
Il985[l and Parker 1797 (|Parkeilll993D . Section [U presents the 
optical spectroscopy of VFTS698, as well as archival opti- 
cal, infrared (IR) and X-ray imaging. The observational 
characteristics of VFTS698 are discussed in Sect. |31 in ad- 
dition to a preliminary analysis. VFTS 698 has previously 



Table 1. Observational spectroscopy for VFTS698. The 
number of epochs associated with each setting is given in 
parentheses. 



been characterised as a B3 la star by Walborn fc Bladei 
(|1997ll and "early Be" bv lBosch et all (Il999ll . but in Sect. 



13.31 we present our own classification based on the proper- 
ties determined here. Section U documents the associated 
atmospheric parameters, including an estimate of the ni- 
trogen abundance. We compare VFTS698 with a number 
of other sgB[e] stars and interacting binary objects in Sect. 
O as well as presenting a plausible evolutionary scenario. 
Our conclusions are drawn up in Sect. |6l 



2. Observations 

2.1. Optical spectroscopy 

Spectroscopy of VFTS698 was obtained as part of the 
Tarantula Survey, using the Fibre Large Array Multi- 
Element Spectrograph (FLAMES) on the Very Large 
Telescope (VLT) with the Giraffe spectrograph. The star 
was observed as part of 'Field B' of the survey, for which the 
observational epochs are listed in the Appendix of Paper I. 
Three of Giraffe wavelength settings were used, as sum- 
marised in Table [T] 

The data reduction was discussed at length in Paper I. 
In brief, the data were reduced using the European 
Southern Observatory (ESO) Common Pipeline Library 
FLAMES package, and consisted of flat-fielding, bias sub- 
traction and wavelength calibration. Separate fibres were 
used to observe the sky in each exposure, with the aver- 
age sky spectrum then subtracted from each of the science 
targets. However, as discussed in Paper I, this often leads 
to imperfect nebular subtraction. In particular, VFTS698 
is embedded in a H ii region with strong nebular emission 
lines which vary over small spatial scales, so accurate neb- 
ular subtraction was not possible. 

Cross-talk between fibres on the array can be a source 
of contamination between objects (see Paper I). This was 
investigated for VFTS698 using the LR03 frames. An in- 
crease of less than 1% in the typical inter- fibre flux was 
found between VFTS698 and one of the adjacent stars on 
the array (VFTS591). Such a small contribution to the flux 
of VFTS698 should not signiflcantly affect our analysis. 

The extracted spectra were normalised over the entire 
spectral region for each setting, or, for features of partic- 
ular interest, in smaller wavelength regions. On occasions 
the richness of the spectra make it difficult to identify line- 
free regions - the implications of this for the silicon spec- 
tra are discussed in Sect. 01 Individual exposures were nor- 
malised using low-order polynomials, with a sigma-clipping 
algorithm used to exclude cosmic-rays within defined con- 
tinuum windows. Spectra within an epoch and, where ap- 



Giraflte set-up 


A- range 


R 


No. Exposures 




(A) 






LR02 


3960 - 4560 


7,000 


15 (5) 


LR03 


4505 - 5050 


8,500 


6 (1) 


HR15N 


6470 - 6790 


16,000 


4(1) 



propriate, between epochs were combined using a weighted 
average together with sigma-clipping. Adopting a median 
spectrum leads to very similar results indicating that the 
choice of algorithm for combining individual exposures was 
unlikely to be a significant source of error. 



2.2. Ptiotometry 

Photometry of VFTS698 is available from the Optical 
Gravitational Lensin g Experiment III (OGLE-III) survey 
(jUdalski et al.l l2008l and references therein), comprising 
multi-epoch observations over eight seasons from October 
2001 to Aprfi 2009. /-band photometry for VFTS698 was 
obtained from the OGLE database (see lUdalskil l2003l for 
further details) and is shown in Fig. [1] in which the mean 
and median measurements for each observing season are 
over-plotted. Note that most of the FLAMES spectroscopy 
is contemporaneous with the last block of OGLE-III data. 

We have supplemented the OGLE-III photometry with 
y-band observations from the 2-m Faulkes Telescope South 
(FTS), obtained as part of a monitoring programme of 
30 Dor during 2009 and 2010. The observations employed 
the FTS Merope camera, which has a field-of-view of 
4.'7 X 4'7 and a (rebinned) pixel size of O'.'278/pixel. Images 
obtained in seeing of greater than 2" were rejected, leaving 
39 observations of the field including VFTS698. Differential 
aperture photometry was obtained for VFTS698 via com- 
parison with five (non- variable) stars of comparable bright- 
ness in the same field. 

To investigate its near-IR behaviour, the Spitzer 
Archiv43 was searched for photometry of VFTS698 and 
nearby sources with comparable y-band magnitudes. 
Observations with the 3.6, 4.5, 5.8 and 8.0/zm filters 
are available from t he Spitzer Legacy Science Program 
(jMeixner et al.l l2006l . SAGE), with photon ietry from the 
survey for known massive stars presented bv lBonanos et all 
(200^. 

Four B-type stars with similar y-band magnitudes (in- 
cluding VFTS450, an 09.7 star) are located within 2' of 
VFTS698. While VFTS698 was detected in aU fihers, this 
was not the case for the comparison objects. For common 
detections, VFTS698 was found to be on average 2 mag 
brighter at 3.6/im, 1.5 mag at 4.5/im, 2 mag at 5.8/im, and 
the only object detected at 8.0^m. Representative visual 
and IR magnitudes are summarised in Table [51 although 
we emphasise that VFTS698 is photometrically variable. 



http:/ /irsa.ipac. caltech.edu/ 
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Fig. 1. /-band light-curve of OGLE-III data from 8 years 
of observations. VFTS698 shows a short-term variation of 
~0.2 mag, increasing up to ^0.5 mag with time. Sohd circles 
and squares show the mean and median values, respectively, 
for each block. The typical uncertainty associated with the 
OGLE-III data is ±0.015 mag, as illustrated in the upper 
left of the figure. 



2.3. X-ray data 

We searched the Chandra Data Archival for images of 
the region including VFTS698. Archival images with the 
ACIS-I instrument of this part of 30 Dor are available 
from 21 Sept 1999, and 21, 22 and 30 Sept 2006. The pri- 
mary data products from the Chandra X-ray Centre's au- 
tomated processing were downloaded and examined, but no 
obvious source was found coincident with VFTS698. This 
supports the analysis of the 1999 data by iTownslev et al.l 
(|2006D who did not report a detection at the position of 
VFTS698 (given a lower detection limit of 5 counts in 
21870 s, which equates to a limiting full-band luminosity 
of ~1 X 10'^'^ ergs s~^ , with the brightest source detected be- 
ing KlO'^^ergss"^). 

We reprojected all of the available individual exposures 
to a common pointing and stacked them to create a deep 
image with an effective exposure time of 120 ks. The deep, 
merged image (Fig. [2]) was filtered to only include pho- 
tons with energies in the range of 0.3 to 10 keV and a 
wavelet algorithm was used to detect point sources. Even 
when an extremely low significance threshold was used (cor- 
responding to a false positive rate among the detections of 
20%), no source was found to be coincident with, or close 
to, VFTS698. 

3. Observational properties and preliminary analysis 

3.1. Optical photometry 

From Fig. [T] there is clear evidence for photometric variabil- 
ity of ~0.6mag on timescales ranging from days to years. 
Two separate trends can be identified: short-term variabil- 
ity with a period of 12.7 ± 0.1 d, and a longer-term variation 
with timescales of at least several hundred days. 




http:/ /cda. harvard.edu/chaser/ 



Fig. 2. Merged Chandra image (smoothed for legibility) 
of the region surrounding VFTS698 (marked by the white 
crosshairs) . 



The short-term variation is illustrated in Fig. [31 Within 
the observational uncertainty the period appears to be con- 
stant between epochs, but the amplitude of the periodic- 
ity increases in magnitude over the 8 years, from ~0.10 to 
'-^O. 25 mag. Each observation block was analysed individu- 
ally by means of a Lomb-Scargie periodogram, confirming 
the consistency of the 12.7 d period. 

The long-term trend appears to consist of a steady in- 
crease in the mean J-band magnitude of ~0.3mag over a 
period of some years, which may have reached a maximum 
at MJD ^54500 (see Fig. [T]). To investigate any possible 
periodicity in this trend, a Lomb-Scargle periodogram was 
produced from the /-band photometry, but now restricted 
to a period range of 100 to 1000 d. The result of this indi- 
cated a period of 400 ± 21 d, but upon folding the period 
to the complete /-band observations, no periodic phase is 
apparent (left-hand panel of Fig. |4]). The difficulty in iden- 
tifying any periodicity in this long-term variation is due 
to the large brightness changes associated with the 12.7 d 
period. When the mean magnitudes for each block are plot- 
ted against the median phase a more convincing variation 
is found, as shown by the right-hand panel of Fig. 2] 

As discussed in Sect. 12.21 VFTS698 was also observed 
over an 11 month period in the V"-band with the FTS. These 
data were found to be in good agreement with the OGLE- 
III data; a Lomb-Scargle periodogram yielded a period of 
12.7 ± 0.1 d (hereafter referred to as the 12.7 d period) with 
a photometric amplitude of ^0.25 mag. The phased light- 
curve from the FTS data is shown in Fig. [5] 

3.2. Infrared photometry 

The near- and mid-IR colours of VFTS698 are compared 
to other massive stars in the LMC in Figs. [5] and [71 
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Fig. 3. OGLE-III /-band photometry of VFTS698 folded with a 12.7 d period. Left-hand panel: First four blocks of 
observations; Right-hand panel: blocks five to eight. For clarity two periods have been plotted for each observation block. 



near-I R photometry is f rom the InfraRed Survey Facility 
(IRSF, i Kato et al.l[2007h and longer wavelengths are from 
the Spitzer SAGE survey. To investigate its behaviour, we 
also include results for VFTS1003, the B[e] star discussed 
in Paper I. As shown in the figures, both VFTS698 and 
VFTS1003 display significant IR excesses over 'normal' 
OBA-type supergiants. 

VFTS698 is clearly displaced from the region of the 
colour/colour plots occupied by sgB[e] stars. At mid-IR 
wavelengths it is co-located with other massive stars such as 
LBVs and Wolf-Rayet (WR) stars, which are known to sup- 
port strong dense winds resulting in copious thermal free- 
free emission at near- to mid-IR wavelengthfl The near-IR 
colours of VFTS698 are bluer than bona fide sgB[e] stars, 
but they indicate a greater continuum excess than found for 
either the LBVs or WRs. The origin of the IR continuum 
emission in VFTS698 is therefore uncertain; in Sect. [S] we 
argue that it appears likely to arise (at least partially) in 
an ionised stellar wind or gaseous circumstellar/binary disc, 
although an additional component from hot dust would also 
appear possible. 



* We note that the source of such excesses in the near-IR for 
LBVs may also be s mall ejecta nebulae surrounding the star 
ijBonanos et al.ll2009f l. 



In contrast, we note that VFTS1003 is co-located with 
the sgB[e] stars in Figs. [5] and [71 suggesting the presence 
of a warm (:~600K, see Sect. [5]) circumstellar dust. Note 
that there is some ambiguity regarding the classification of 
VFTS1003 in Paper I, with it potentially a sub-luminous 
sgB[e] or an over luminous, pre- main sequence 'Herbig' B[e] 
star. The near- and mid-IR photometry would appear to 
support the sgB [e] classification, but the close proximity of 
VFTS1003 to the central region of 30-Dor casts doubt on 
the reliability of the Spitzer photometry. 

3.3. Spectral classification and luminosity 

VFTS698 shows a rich absorption- an d emission-lin e spec- 
trum, as illustrated in Figs. [8l and ITOl Tamers et al.l (|l998f) 
presented an overview of the sub-classification of B[e] stars, 
from which we classify VFTS698 as a B[e] supergiant. Our 
two primary criteria were the presence of forbidden emis- 
sion lines and a large luminosity. In addition, five secondary 
criteria were considered: indications of mass-loss from 
double-peaked Balmer emission lines; presence of narrow, 
low-excitation emission lines and broad, higher-excitation 
absorption lines; a strong nitrogen spectrum; a large extinc- 
tion; photometric variability. For the last point, VFTS698 
shows larger variations than most sgB[e] objects (typically 
0.1 to 0.2 mag), but it is not unique - e.g. R4 in the SMC 
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0.0 0.2 0.4 0.6 0.8 1.0 



Fig. 4. Left-hand panel: /-band OGLE-III photometry of VFTS698 folded to a 400 d period. Right-hand panel: The 
median magnitude and date for each block of OGLE-III /-band data as a function of phase. The bar indicates the typical 
phase coverage by each block. 



Table 2. Photometry of VFTS698 and four nearby, early B-type (or late O) VFTS sources with comparable V"-band 
magnitudes. Photometry of VFTS1003, provisionally classified as B[e] in Paper I, is also included for comparison. 



Star 


Alias 


Lit. 


VFTS 


V 


B-V 


J 


H 


Ks 


[3.61 


[4.51 


[5.8J 


[8.0] 






Sp. type 


Sp. type 




















VFTS4501' 


Mk50 


0N9: I' 


09.7 


13.60 


0.20 


13.08 


12.91 


12.89 


11.00 


10.65 


10.52 




VFTS6521' 


Mk5 


B2 Ib:^ 




13.88 


0.20 


13.40 


13.28 


13.22 


12.97 




12.73 




VFTS698* 


Mk2 


B3 la^ 




13.68 


0.44 


12.34 


11.89 


11.43 


10.60 


10.23 


10.02 


9.66 


VFTS732* 


Mkl 






13.03 


0.19 


12.35 


12.19 


12.08 


11.85 


11.73 


11.30 




VFTS733* 


P1988 


B0.5 




14.28 


0.12 


13.90 


13.82 


13.74 


13.59 








VFTSlOOSt 


S99-283 




Bfel? 


16.10 


0.23 


13.86 


13.32 


12.45 


10.49 


9.73 


9.02 


8.34 



N otes. Spectral types from: IWalborn fc BladesI (|l997l h^^' iBosch et all (|1999I 1. Optical phot ometry from: IParked (|1993D : 
ISelman et aD (|l999l '). Near-IR ( JHKs) pho tometry from IRSF Magellanic Clouds catalogue (|Kato et al.ll2007f ). Mid-IR pho- 
tometry taken from Spitzer SAGE ijMeixner et al...2006. ). 



displays comparable variability (jZickgraf et al.l [19961 ). The 
strong photometric variability and mid-IR colours appear 
to weaken the sgB[e] classification, but there remains suffi- 
cient evidence to justify the classification as an early-type 
supergiant with '[e]' phenomena. 

The photometric variability of VFTS698 makes its lumi- 
nosity difficult to estimate. However we can use a number of 
constraints to obtain a best estimate. Using the photometry 
from Table[2j ad opting ( B—V )o = —0.24 mag and a distance 
modulus of 18.5 (|Alvesll2004[) to the LMC, we obtain E{B- 
y) = 0.6mag and My = —6. 7 mag. Adopting an appropri- 
ate bolometric correction for an ea rly B-type supergiant 
(—2.4 mag, see lCrowther et al]|2006l ) results in a bolomet- 
ric magnitude of —9.1 mag, equivalent to log L/Lq = 5.6. 

3.4. Spectral characteristics 

We now discuss the FLAMES spectra, which appear to have 
three distinct components: 

— Characteristic sgB[e] spectrum: Low-excitation permit- 
ted and forbidden lines from iron-group elements; 



— Cool component: Absorption lines from species consis- 
tent with a mid-B spectral classification, e.g., Hei and 
Sin; 

— Hot component: High-ionisation stages (Niii, Siiv and 
Hen), consistent with an early B- or late 0-type spec- 
trum. 

The apparently cool features display small RV varia- 
tions, with far larger variations in the lines attributable to 
a hot component. The RV variations of example metal and 
helium lines (Siiv, Niii and Hel) are shown in Fig. [9l We 
now consider these different components in turn, although 
there is inevitably some overlap in the discussion. 

3.4.1. Characteristic sgB[e] spectrum 

To inform the description of VFTS698, Fig. [TU] compares 
the LR02 data with spectra of two sgB[e] stars in the 
SMC, obtained with the Cassegrain Echelle Spectrograph 
(CASPEC) on the ESQ 3.6-m telescope (degraded and re- 
binned to match the resolution of FLAMES). Zickgr af et al.l 
()1996f ) noted the prominent emission of the Balmer and Fe ii 
lines (both permitted and forbidden) in R4. Similar Fe ii 
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_ ♦ LMC sgB[e] 

» LMC LBV+WB9-1 1 h 

X LMC WN (excL WN9-1 1 h) 

a LMC WC 

o LMC DBA supergiant 



K,-[8.0] 

Fig. 6. J - Ks colour vs. Ks-[S.Q\ colour for VFTS698, 
other emissfon-line objects and OBA-type supergiants in 
the LMC. The recently discovered B[e] star VFTS1003 is 
also shown for comparison. Note that strong emission lines 
are also expected to contribute to the near-IR continuum 
emission of WC stars. 



-0,5 O.C 



1,0 1.: 



Fig. 5. Upper panel: Differential F-band photometry of 
VFTS698 from 11 months of FTS observations (open cir- 
cles) . Mean and median magnitudes for the two seasons are 
indicated by solid circles and squares, respectively. Lower 
panel: FTS photometry folded to a 12.7 d period. 



features are seen in VFTS698, but they are accompanied 
by narrow absorption lines (see Fig. [TT|) : such structure is 
also seen in the Tin and Cm lines in VFTS698. This two- 
comp onent structure is n ot seen in R4, but is present in 
R50 (jZickgraf et al.lll98fil upper spectrum in Fig. ITU)) . 

The most prominent permitted lines of Fell are due 
to the multiplets 27, 28, 37 and 38, with an excita- 
tion potential '-■^2/7-5.6eV (similar to those identified by 
[Zickgraf et al.lll996[ ). whilst for Tin and Cm, they are due 
to multiplets 41 and 31, with an excitation potential ~1.2- 
4.0 eV and ~3.9-6.8eV, respectively. The forbidden tran- 
sitions mainly originate from [Fell] multiplets 4F, 6F, 7F, 
20F and 21F, with excitation potentials of 0.1 -3.0 eV. A 
number of other forbidden lines are also detected, including 
[Feiii] at 4658, 4701, and 4733A, [Nii] at 6548 and 6583A, 
and [Sii] at 6717 and 6730A. However, it is unclear if the 
forbidden transitions of [N ii] and [S ii] originate from ma- 
terial directly associated with VFTS698 as these features 
are also observed in the FLAMES spectra of other nearby 
objects, where they originate from the strong nebular emis- 
sion in 30 Dor (Paper I). 

The narrow-lined, iron-group absorption features have 
been fitted with Gaussian profiles to estimate RVs. This 



Fig. 7. Mid-IR colours for VFTS698, other emission-line 
objects and OBA-type supergiants in the LMC. VFTS698 
appears to have an IR-excess greater than normal OBA 
supergiants, but not to the same extent as other sgB[e] 
stars. 



was performed for all epochs of LR02 spectroscopy, and no 
evidence for RV variations was found, with a mean RV of 
259kms~^ and a sample standard deviation of ±15kms~^ . 
This spectral component appears to lie close to the rest 
frame of 30 Dor, as measured from over two hundred B- 
type stars (<RV> = 270±17kms~^ , Kennedy et al. to be 
submitted). The low excitation potentials of the iron-group 
spectrum, together with the RV measurements which in- 
dicate a static environment, suggest they originate within 
low density, circumstellar material. 
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Fig. 8. Optical spectroscopy of the HR15N observations of 
VFTS698. See Fig. [10] for details on the LR02 and LR03 
settings. 



ably originate from the same source. The He i transitions 
show consistent RV measurements within an epoch (with a 
sample standard deviation of ±15kms~^ ) but with a vari- 
ation of 33 ± 11 km s~^ over the five spectroscopic epochs. 
The Si II and Mg ii lines were subject to blending that made 
accurate estimates difficult, but both species were found to 
have RVs within typically ^20 km s~^ of those of the Hei 
lines. If we consider these lines to originate from a sin- 
gle plasma, this could be the B[e] supergiant photosphere 
with a mid-B spectral classification. Additionally a Lomb- 
Scargle periodogram for the He i RV estimates returned a 
periodicity of 12.7 ± 0.1 d, implying a correlation with the 
photometric variability. 

In Fig. [T2] the strongest He I diffuse triplet transitions, 
shifted to their rest wavelength, are shown as a function 
of their phase, adopting the 12.7 d photometric period. It 
should be noted these features are susceptible to contamina- 
tion by narrow nebular emission, which may not have been 
completely removed. However, variations in the P Cygni na- 
ture of their profiles are clearly visible. Hereafter we refer 
to these features as the 'cool stellar' component. 




1130 1135 1110 1115 1085 1090 1095 1100 1105 

Wovelength (A) Wovelerglh (A) 

Fig. 9. Example RV variations displayed by helium and 
metal lines. Fitted Gaussian profiles are over-plotted and 
rest wavelengths are indicated by the vertical dashed lines. 
Left-hand panel: He I 4437A, in which the 'hot', second com- 
ponent is noted by the marks above the continuum. Right- 
hand panel: Si IV 4089A and Niii 4097A. 



3.4.2. Cool component 

Absorption lines from Nil, Sill and Mgii are observed in 
the spectrum of VFTS698 with profiles that are broader 
than those from the iron-group species discussed above. 
There is also a strong neutral helium spectrum, unusual 
for objects exhibiting B[e] features, which show a two- 
component structure (see Fig. [3]). 

As discussed below the Hel diffuse triplet lines show 
varying P Cygni profiles but the singlet lines appear reason- 
ably symmetric. Determinations of the RVs of the stronger 
component of these singlet lines and of the singly-ionised 
silicon and magnesium transitions indicate that they prob- 



3.4.3. Hot component 

High-ionisation absorption lines were also identified in the 
spectrum, viz. Sim, Si IV and Niii features as well the 
strongest He li transition in this wavelength region, 4686A. 
The Si IV and Niii features were observed in the multi- 
epoch LR02 data, which reveal large RV variations with 
a range of ^^260 km s^^. Including the LR03 observation 
(which includes the Hen, and Sim lines) this range in- 
creases to '^420kms~^ . The weaker component of the Hei 
spectrum (discussed in the previous section) has a simi- 
lar velocity structure to these high-ionisation metal lines. 
Additionally, the N ii spectrum from the LR03 observation 
also displays a two component structure, with the RV esti- 
mate of one component compatible with those of the other 
high-ionisation features. Although the weaker He i compo- 
nent appears to have a similar RV variation to those of the 
high-ionisation metal lines, it is normally blended. Hence, 
the velocity range was determined only from the metal lines. 
Assuming that these lines all originate from a single stellar 
object, it would have an early-B/late-0 spectral classifica- 
tion. 

In Fig. [T3] the RV estimates from the high-ionisation 
features are folded with the 12.7 d period found from the 
OGLE-HI data. Also shown are the OGLE-HI data for the 
block that overlapped in time with most of the FLAMES 
spectroscopy. From the figure the RV data appear to be 
consistent with the variations in magnitude and, although 
there are insufficient points to obtain an accurate period, a 
Lomb-Scargle periodogram returns a period of 10.5 ± 0.1 d, 
which appears consistent with the photometry. 

A detailed list of the absorption and emission features 
for VFTS698 are presented in Table [8] (available online). 
Values for the central wavelengths (corrected for the RV 
estimated from the iron-group spectra), equivalent widths, 
and full-width half-maxima of the profiles are provided. For 
spectral lines observed with the LR02 setting, these mea- 
surements are from the first epoch. 
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Fig. 10. The combined LR02 and LR03 spectra of VF TS698 (lower spectrum), compared w ith two sgBfe] stars 
in the SMC: R50 (upper spectrum, from iZickgraf et al.l [l986i ) and R4 (middle spectrum, from IZickgraf et al.l 119961 : 
iPodsiadlowski et alJ[2006h . 



Table 3. Radial velocity measurements of absorption features for selected transitions given for both hot and cool 
components. Measurements quoted in kms~^ . 



Hot Components 



Cool Components 



MJD 


He I 


Hen 


Sim 


Si IV 


Nil 


Niii 


He I 


Sin 


Sim 


Nil 


Mgii 


54804 


-84±14 






-86±11 


-99±20 


-91±4 


-7±10 


4±11 




-40±20 


-14±10 


54808 




267±10 


279±10 




287±10 








54±10 


-4± 3 




54836 








-12±11 




-22±17 


0±15 






15±10 




54867 


-106±24 






-129±10 


-145±20 


-110±26 


-20±12 


-5± 9 




-67±20 


-14±13 


55108 


-96±29 






-142±11 


-178±20 


-136±16 


-13±13 


-15±12 




-66±20 


-29±10 


55114 


103±21 






130±10 


56±20 


100±13 


17±11 


18±11 




7±20 


47±10 



3.5. Projected rotational velocities 



Assuming that the hot and cool components originate in 
stellar photospheres, we can attempt to estimate their pro- 
jected rotational velocities. We hav e used a Fourier trans- 
form ( FT) approach as discussed bv lSimon-Diaz fc Herrerj 
(|2007f) , supplemented by fitting rotationally- broadened the- 
oretical line profiles (PF) to the observed spectral features. 
For the latter we assumed that rotation dominated other 
broadening mechanisms and hence these are best consid- 
ered as upper limits. Further details on the pro file-fitting 
methodology can be found in lRvans et al.l (|2002[ ). 



For the hot component we used the Siiv lines at 4089 
and 4116A and the Niii line at 4097A. For those LR02 
epochs where the lines were discernible from features in 
the cool or circumstellar spectra, the data were shifted to 
a common rest frame and combined. This led to relatively 
clean profiles for the Siiv lines but the Niii line was still 
affected on its red wing by emission. The estimates are 
listed in Table U and are in reasonable agreement for both 
methodologies. Additionally, as expected, the PF estimates 
(which are upper limits) are generally larger than those 
from the FT approach. This provides indirect support for 
the validity of the latter. A simple average of the FT results 
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Fig. 12. The Hci 4026 and 447lA normalised line profiles (left- and right-hand panels, respectively) shifted to their rest 
wavelengths, with their continuum level indicating their phase given a single 12.7 d period. The three epochs at a phase 
of 0.2 display similar structure, whereas later in the period (at phases of 0.6 and 0.8) there are significant changes in 
morphology. 
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Fig. 11. Fell , Cm and Tin emission, and shell- like fea- 
tures from the (first epoch) LR02 spectrum of VFTS698. 



leads to a value of 66 km s"'^ for the hot component, with 
the spread of estimates and the upper limit from the PF 
estimates implying a conservative error bar of ±10kms~^ . 

For the cool component, lines due to several metal 
species were identified. These suffered from blending with 
cither the hot or circumstellar components. This was ex- 
acerbated by the features having a small RV variation 



1 3.5 -800 

-0.5 0.0 0.5 1.0 1.5 2.0 2.5 

Phase 

Fig. 13. Upper panel: Radial velocities of the hot and cool 
components (open and solid circles, respectively) folded 
with a 12.7 d period. Lower panel: Single block of phased 
OGLE-III data. The radial velocity measurements for the 
hot component appear to be consistent with the periodicity 
of the photometric data. 



and therefore not shifting with respect to the circumstel- 
lar features between epochs. Hence the stronger Hei fea- 
tures were used, although their analysis is complicated 
by their greater intrinsic broadening. Additionally, several 
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Table 4. Estimates of the projected rotational velocity 
(usini) of the cool and hot components from the profile 
fitting (PF) and Fourier transform (FT) methods. 



Component 


Feature 


PF 
(kms~ 


vsini 

FT 

1) (kms-i) 


Hot 


Si IV 4089A 


79 


63 




Si IV 4116A 


68 


61 




N III 4097A 


77 


74 


Cool 


He I 4009A 


70 


60 




He I 4143A 


64 


50 




He I 4437A 


58 


49 



lines showed P Cygni profiles (see Sect. 13.4. 2p and were 
not considered. The estimates are summarised in Table 2] 
They exhibit the expected systematic shift between the PF 
and FT approaches and lead to a mean FT estimate of 
53±10kms~^ for the cool component; the error estimate 
again arises from the spread of values and the upper limit 
implied by the PF results. 

4. Atmospheric parameters 

The complexity of the VFTS698 spectrum, containing sev- 
eral different components, makes any quantitative analysis 
both difficult and prone to systematic error. However the re- 
sulting quantitative information could be useful despite the 
large uncertainties. We have therefore attempted to obtain 
constraints on the effective temperature, surface gravity 
and nitrogen abundances for the two comp onents using th e 
non-LTE model atmosphere grid tlusty (Hubenv 1988'); 
furthe r det ails of the gr i ds ca n be found in iRvans et al.i 
(pool and lDufton etall (l2005l) . 

4.1. Hot component 

The well-resolved spectral features of the hot stellar com- 
ponent made estimating its effective temperature from the 
Si Iii/Siiv ionisation equilibrium relatively straightforward. 
Due to the strong Fe 11 emission and absorption spectrum 
between 4500 and 4560A (shown in Fig. ITO)) , as well as the 
wavelength cut-off at 4560A, the Sim lines were not mea- 
surable in the LR02 setting. The longer wavelength cover- 
age of the LR03 setting and the radial velocity of the hot 
component (>200kms~'^ ) in these data allowed equivalent- 
width measurements of Sim 4567Afor features that had a 
RV consistent with the hot and cool photospheric spec- 
tra. The equivalent widths of Siiv were taken from mea- 
surements of the 4089 and 41 16A lines, in the LR02 data. 
The two Si iv lines displayed variations in intensity between 
epochs of 12% and 26%, respectively. The average equiva- 
lent widths for both Siiv lines gave an effective tempera- 
ture estimate of ~27000K. However, the equivalent widths 
will be underestimated due to contributions from the other 
components in the spectrum. We assumed that this contri- 
bution was the same for all lines, which is not unreasonable 
given the relatively small wavelength separation. Then we 
investigated its effect on the effective temperature estimate 



by scaling the equivalent widths. Doubling the equivalent 
widths led to an increase in the estimate of 1 000 K. 

A lower limit to the effective temperature of the hot 
component was estimated from the Hen absorption line 
at 4686A. Rotationally-broadened profiles with a f sin? of 
70 km s~^ were generated from the tlusty grid of model 
atmospheres. From a fit an effective temperature of ~ 
24 000 K was found together with a gravity, log 3.3. If 
the hot spectra contributed 50% of the continuum flux this 
estimate would be increased by 1 500 K. The uncertainty in 
these estimates is difficult to judge. A lower limit of 24 OOOK 
would appear reasonable both from the He 11 spectrum and 
the presence of strong Si iv lines. An upper limit of 28 000 K 
is reasonable from the well observed Si ill. Hence we adopt 
26 000 ± 2 OOOK, but acknowledge that there remains the 
possibility of additional systematic errors from our assump- 
tion of a classical non-LTE photospheric model. 



4.2. Cool component 

For the cool photospheric spectrum, the detection of Sill 
and Si III lines also allows the silicon ionisation balance to 
be used to estimate the effective temperature. Equivalent- 
width estimates were obtained for the Sill 4128 and 4130A 
absorption lines from the LR02 data, and for Sim at 
4567A. Those for the former were difficult to estimate as 
the RVs were inconsistent within an epoch (see Table [3]). 
This discrepancy is probably due to a weak emission fea- 
ture around 4129A, and a Cm absorption component seen 
to be blended with the Sin 4130A line. Assuming there is 
no flux contamination for the other spectral components, 
this leads to an effective temperature estimate of 17 800 K, 
for an assumed logarithmic gravity of log q= 3.0 based o n 
other B-type supergiants taken from iHunter et all (|2008[ ). 
Making similar assumptions as in the analysis of the hot 
component, i.e. doubling the equivalent widths, leads to 
an estimate of 18 400 K. Again the uncertainty is difficult 
to estimate. However, the clear presence of both ionisation 
states would imply that, excluding unknown systematic er- 
rors, an error of ±2 000 K would appear to be appropriate 
and therefore we adopt a Toff of 18 000 ± 2 OOOK. 

Apart from a weak constraint from the profile fitting of 
the Hen line, it was not possible to estimate the surface 
gravity of either components due to the strong emission in 
the B aimer lines. 

In Sect. 13.31 a limit was placed on the luminosity of 
VFTS698 assuming it to be a single source. However, in- 
formed by the estimates of effective temperature further 
constraints could be placed on each of the components. 
Both the hot and cool spectra appear to be prominent, 
hence their contribution to the overall spectrum is thought 
to be comparable. Based on the effective temperature esti- 
mates from Table [B] we deduce luminosities for the hot and 
cool components of log L/Lq = 5.3 and 5.0, respectively , 
using bolometric corrections from (jCrowther et al.l l2006f) 
and identical absolute visual magnitudes. To consider these 
components as main sequence objects would lead to a mass 
ratio estimate of ^2, but given our analysis of the VFTS698 
system classifies it as a sgB[e] system this consideration is 
not applicable. 
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4.3. Nitrogen abundance estimates 

The photospheric abundance of nitrogen is important for 
constraining the amount of nucleosynthe ticany-processed 
material present in a stellar atmosphere (|Heger fc Langeil 
l2000HPrzvbilla et al.ll2010HBrott et al.ll201l[ ) As Nii lines 
from both stellar components could be identified in the 
spectra, we can set lower limits on their nitrogen abun- 
dances. 

For the hot component, adopting an effective tempera- 
ture of 26 000 K, a gravity of log g= 3.3 and a microtur- 
bulence of lOkms"^ leads to an estimate of ^7.3 dex from 
three N li lines at 4600A. Note that adopting quantities for 
the gravity and microturbulence that would be expected for 
a star on the main sequence results in an effective tempera- 
ture estimate of 30 000 K and a surface nitrogen abundance 
of '^S.Odex. These weaker transitions were initially chosen 
because of their large velocity separation. The stronger N ii 
line at 3995A was identified in all the LR02 epochs, but 
the blending of the two stellar signatures made it difficult 
to disentangle the equivalent width estimates. Following a 
similar procedure to that used in measuring the RVs of 
the He I profiles, two-component Gaussian fits were em- 
ployed, with the equivalent width of the hot component 
taken to be that with a RV component corresponding to the 
higher-excitation metal lines. Although there was a varia- 
tion of the line strength between epochs, a mean value of 
50 mA implied a nitrogen abundance of 7.2 dex using the 
same atmospheric parameters. 

The above analysis assumes that all the continuum fiux 
arises from the hot photospheric spectrum and therefore 
abundances were also estimated assuming 50 and 25% of 
the flux came from the hot component, as summarised in 
Table [SJ The fact that both spectra are clearly visible would 
suggest each component might contribute ^^50% to the to- 
tal flux, which would then suggest a nitrogen abundance 
of 8.0 dex for the hot component. However, regardless of 
the fiux contribution, there would appear to be an enrich- 
ment of nitrogen as, at the limit of a 100% contribution, 
the estimates show a 0.6 dex enhancement. 

The same methodology has been used to estimate a ni- 
trogen abundance for the cool spectral features although in 
some cases the equivalent-width measurements lay outside 
the theoretical grid at the 25% fiux level. The equivalent- 
width estimates of the 3995A line were found to vary more 
between epochs than for the hot spectra but, if we adopt 
the average value, we estimate the lower limit to the ni- 
trogen abundance as 7.9 dex, i.e. an enrichment of 1 dex. 
Adopting a 50% fiux contribution leads to an estimate of 
~8.6dex. 

Adopting a 50% flux contribution from the hot and cool 
components leads to silicon abundances between 7.0 and 7.2 
dex which are in good agreement with those found for early- 
type stars in the LMC (see, e.g.. lHunter et al.ll2008[) . pro- 
viding some indirect support for the methodology adopted. 
However we stress that the estimates of the atmospheric pa- 
rameters and nitrogen abundances should be treated with 
caution because of the complexity of the VFTS698 spectra. 
In order to discuss the evolutionary status of VFTS698 in 
the following sections, we provide a summary table of the 
estimated properties of both hot and cool spectra in Table 

m 



5. Discussion 

The rich emission-line spectrum of VFTS698, dominated by 
Hi, Fell and [Feii], is characteristic of stars demonstrating 
the B[e] phenomenon. Lanicrs ct al. (1998) describe this be- 
haviour as being present at a number of evolutionary stages 
in the life-cycle of both high- and low-mass single stars and 
binaries; essentially corresponding to the presence of a rich, 
density-stratified dusty and gaseous circumstellar environ- 
ment, excited by a hot ionising source. We now discuss the 
nature of the circumstellar environment and that of the cen- 
tral 'engine' of VFTS698 from the observational properties 
described in Sect. [3] 

5.1. Circumstellar material 

The differing line profiles of both the Fe ii and [Fe ii] tran- 
sitions (e.g.. Fig. [TUl) imply that they arise in kinematically 
distinct, gaseous regions of the circumstellar environment. 
Indeed, by analogy to classical Be stars, the shell profiles of 
the Fe ii transitions - characterised by deep central absorp- 
tion troughs - suggest absorption in a flattened circumstel- 
lar disc or torus seen close to edge on. This is also analogous 
to the P Cygn i profiles that arise in a spherically expanding 
envelope (e.g. iDachs et al.lll992ll . Following from this, the 
bona fide P Cygni profiles observed in e.g. the He i 4026A 
transition suggest a spherical outflow from the system. Such 
hybrid behaviour is common amongst sgB[e] stars, where 
it is often interpreted as resulting from a 'composite' wind 
consisting of a dense, outflowing circumstellar disc and a 
high-velocity polar wind (e.g. Zickgraf ct al. 1985). 

However, as demonstrated clearly by Figs. [5] and [3 the 
IR propertie s of VFTS698 differ from those of bona fide 
sgB[e] stars. iKastner et al] ()2006l . I2O1O0 showed that the 
IR emission from these stars is consistent with a population 
of warm grains located within a circumstellar disc with an 
inner gap resulting in a deficiency of hot dust. Upon consid- 
eration of the dust chemistry they further speculated that 
these are long-lived quasi-Keplerian discs, with material de- 
posited from a previous red supergiant phase. By analogy to 
similar structures associated with post-AGB binary stars, 
it is possible that these discs are in fact circumbinary with 
an (unseen) companion responsible for the disc truncation. 
In contrast, VFTS698 does not appear to support large 
quantities of cool dust, although the presence of an appar- 
ent near-IR continuum excess raises the possibility of the 
presence of hot dust; detailed modelling will be required to 
characterise the nature of this excess, given that we would 
also expect a significant free-free emission from the gaseous 
component of the circumstellar environment. 

5.2. Photospheric sources 

We next turn , t o the central, exciting star(s) of VFTS698. 
iLamers et al.l (|1998[ ) summarised the stellar properties 
of a representative sample o f Ma gellanic Cloud sgB[e] 
stars, from investigations by IZick graf et al. (1986) and 
iGummersbach et al.l (jl995t) . The temperatures of both our 
putative hot and cool components (Table IH) are consis- 
tent with these values, although the luminosities we infer 
are, for the most part, towards the lower range demon- 
strated (cf. Sect. |4]), with 5.3 and 5.0 dex for the hot and 
cool components, respectively. The luminosity and effective 
temperature estimates of VFTS698 and other B [e] and B- 
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Table 5. Equivalent width (EW) measurements (quoted in mA) and abundance estimates for both Nil components. 



Epoch Flux 


4601 


4607 


4613 


4643 


3995 


(MJD) % 


EW 


N 


EW 


N 


EW 


N 


EW N 


EW N 


Hot component: 


54808 100 


34 


7.6 


23 


7.5 






51 7.8 


50 7.2 


50 


68 


8.0 


46 


7.9 






102 8.3 


100 7.6 


25 


136 


8.6 


92 


8.4 






204 


200 8.3 


Cool component: 


54808 100 


110 


8.3 


120 


8.5 


68 


8.2 




170 7.9 


75 


147 


8.5 


160 


8.6 


91 


8.4 




230 8.2 


50 


220 


8.6 


240 


8.5 


136 


8.6 




340 8.6 



Notes. Each EW was measured assuming different amounts of photospheric flux. For both components the adopted atmospheric 
parameters were taken from Table [B] denotes the EW is outside of the theoretical grid. 



Table 6. Summary of derived parameters for the VFTS698 system, assuming 50% of the flux comes from each component. 



Component 


RV amplitude 


V sin i 


log L/Lq 




log g 




N 


Si 




(kms-^) 


(kms-i) 




(K) 


log (cms~^) 


(kms-^) 






Hot 


210 


66 


5.3 


26 000 


3.3 


10* 


8.0 


7.0 


Cool 


16 


53 


5.0 


18 000 


3.0* 


20* 


8.6 


7.2 



Notes. *■*' Assumed value of log 3 or microturbulance 



3 L 
4,6 



log (Teff) 



Fig. 14. Hertzsprung-Russell diagram comparing the two 
components of VFTS698 ( solid squares) with B -type super- 
giants (open circles) fro m iHunter et all (I2008D and sgB[e] 
stars (solid circles) from lLamers et al.~()1998[ ). 



type supergiants are illustrated on the Hertzsprung-Russell 
(H-R) diagram shown in Fig. [M) For the majority of B [e] 
objects it is evident that both components of VFTS698 ap- 
pear at the lower luminosity end, but are brighter than the 
majority of B-type supergiants. The evolutionary nature 
of these lower luminosity (and mas s) B[e] stars is uncer- 
tain (e.g. iGummersbach et al.lll995h but we note that this 
phenomenon is not uncommon or confined solely to the 
Mage llanic Clouds, with a population also identified within 
M33 (IClark et al.ll2012D . 



A key observational constraint placed on the stellar 
component (s) of VFTS698 is the presence of pronounced 
spectroscopic and photometric variability. The presence 
of variability of th e Siiv lines echoes the findings of 
IZickgraf et all (|1986[ ). who detected Siiv 4089 and 4116A 
absorption in th e spectrum of Hen S22 where previously 
Muratoriol (|1978[ ) observed no such features (although the 
source of this variation remains unclear). More telling is the 
presence of RV shifts in lines associated with both hot and 
cool components. Such behaviour is characteristic of both 
stellar pulsations and binary motion, although the magni- 
tude of the RV changes in the hot component is significantly 
greater tha n expected from pulsations alone (jClark et al.l 
I2OIOL l2012l ). Indeed, with two plasmas with significantly 
different temperatures (supported by the wide range of ion- 
isation states observed, e.g., Siii/Siiii/Siiv and N ii/N ill), 
with very distinct RV behaviours, it is highly plausible that 
VFTS698 is a hierarchical stellar system. 

The multi-year photometric dataset strongly supports 
such an assertion, revealing a persistent 12.7 d periodicity 
superimposed over a long-term variation for a duration of 8 
years fSect. [3T|) . While pulsationally-induced photometric 
variability appears ub iquitous across th e upper reaches of 
the H-R diagram fe.g. lClark et al.ll2012f ) the pulsations ob- 
served in the so-called supergiant a Cygni variables do not 
replicate the strict periodicity observed in VFTS698 over 
the 8 year observational period. Specifically, we may con- 
trast this behaviour to that of the Galactic B[e] star GG 
Car, which Gosset et al. (1984) showed to have a photomet- 
ric period of 31.02 d, with an amplitude of ^0.5 mag that is 
stable o ver a 4 year period. Subsequent spectroscopic obser- 
vations (jGosset et al.lll985l) revealed correlated RV changes 
suggesting a binary nature for this system. Unfortunately, 
with only six epochs of LR02 spectroscopy for VFTS698, we 
are unable to determine accurate periodicities, although, for 
completeness, we note that the RV variations were found 
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in Sect. 13.41 to be reasonably compatible with the 12. 7d 
photometric period. 

5.3. Interacting binary 

Under the assumption of binarity, what may we conclude 
about the properties of the VFTS698 system? Its photomet- 
ric properties sugge st obvious similari ties with the class of 
DPVs identified by iMennickent et all (2003) from OGLE- 
II data. These intrinsically blue objects demonstrate both 
short- ( Pi , days) and long- term periodicities (P2, hundreds 
of days) - IMennickent et al.l suggested that these are related 
to the orbital period and precession of an elliptical accretion 
disc and/or episodic mass loss, respectively, within a semi- 
detached binary systerr0- As such they would be closely 
related in evolutionary terms to the better doc umented 
_W_Serpentis//3 Lyrae lainaries in the Galaxy (e.g. [Tarasovl 
I2OOOI and references therein) . In these short-periocQ, inter- 
acting systems the mass loss from the initially more mas- 
sive component results in the formation of an accretion 
disc/torus (which veils the mass gainer) and in the rever- 
sal of the mass ratio of the system. A consequence of this 
configuration is the presence of photometric modulation on 
the time-scale of the orbital period, as well as binary RV 
variability in the spectrum of the primary. Binary motion 
is also visible in the spectroscopic signature of the accre- 
tion disc around the mass gainer m the subset of systems 
where it may be isolated (cf. lTarasovll2000l ). while some sys- 
tems also demo nstrate long t erm, aperiodic, photometric 
variability (e.g. IStruoatl Il987f). Finally, the mass transfer 
drives a complex circumstellar environme nt, with evidence 
for b oth polar mass-loss (e.g. /3 Lyrae; iHarmanec et al.l 
I1996D and the pres ence of a circumbinary disc (e.g. BY 
Cru; lTarasov ll2000t) in some systems. While the presence 
of P Gygni profiles in selected He i transitions of VFTS698 
is suggestive of the former, the lack of RV variability in the 
Fell/ [Fell] emission-line spectrum argues for dynamical de- 
coupling of the material from which it originates. 

Since VFTS698 appears to demonstrate the majority 
of the physical phenomena associated with W Serpentis- 
type, semi-detached eclipsing binaries, we consider such 
an identification as highly promising. A particularly in- 
teresting comparator is the massive (~7Mq-|-3OM0) bi- 
nary RY Scuti, due to the similar orbital period (~lld 
modulation present in both photometry and spectroscopy) 
and spectral t ype of the mass donor/secondary (09.7 
Ibpe; IWalborn I [1982') compared to the hot component in 
VFTS698 and presence of circumbinary ejecta. Although 
the primary /mass gainer in RY Scuti is enshrouded by 
circumstellar material, u sing tomographic reconstruction 
iGrundstrom et all (j200 7f) estimated a BO. 5 I classification 
for the "massive companion" but suggest that the latter 
is actually the spectrum of the photosphere due to an ac- 
cretion torus directly associated with it. The RV estimates 
for the hot component of RY Scuti were found using high- 
ionisation metal lines (Siiv and Niii), whereas the wings 



^ IMennickent et aP relate the periodicities in DPVs by 
Pi = 3 5.2 ± 0.8 P2; later revised by IMennickent fc Kolaczkowskil 
l|2009n to ~33. From our period search of VFTS698 (Sect. |3TT|| 
we find a 12.7 d period, and some evidence for a 400 d period, 
leading to a ratio of ~31. 

^ E.g.: /? Ly rae (12.9 d), W Serpentis (14.16 d), V367 Cyg 
(18.16 d), from lTaraso^ (|2000l l. 



of Ha, He II 4686 A emission, and Si ill absorption were the 
only features capable of reproducing the anti-phase for the 
cool component. Unfortunately, due to the singular LR03 
and HR15N observations for VFTS698, an equivalent RV 
analysis is not possible. Furthermore, as discussed in Sect, 
m Si III was detected for the cool component but because of 
rich Fe 11 emission around 4550A, multi-epoch data of this 
ionisation state was unavailable. Nevertheless, by analogy 
to RY Scuti it might be supposed that the spectroscopic 
features of the cool companion in VFTS698 actually arise 
in the photosphere of an accretion disc/torus around the 
veiled mass gainer. In such a scenario one would anticipate 
that RV changes be anti-correlated with those observed in 
the hot component, an effect that does not appear to be 
present in our limited dataset. However, we note that if 
the spectrum of the cool component was formed at the LI 
Lagrange point (where Roche lobe overflow will occur - pre- 
sumably as the material will be at its most dense there), 
then this would lie inside the centre of mass of the system, 
and so we would expect to see a weak correlation between 
the hot and cool component spectra, as is observed. 

Further RV observations of VFTS698 are required to in- 
vestigate this discrepancy. If confirmed, one possible expla- 
nation might be the presence of a third body in the system, 
although in this case one would expect a significantly longer 
period than the provisional ^12.7±0.1d period determined 
from the current dataset. Moreover, in comparison to RY 
Scuti, the richer Fe 11/ [Fe 11] emission-line spectrum suggests 
differences in the geometry/composition of the circumstel- 
lar/binary envelope. 

Although no B[e] phenomena has been detected in the 
comparison interacting binariefl it's suggested the SMC 
sgB[e] R4 is a binary merger, the product of a cl ose bi- 
nary interaction (Zickgraf et al. 1996; Podsiadlowsk i et al.l 
120061) . iLanger fc Hegeii (|1998il have attempted to explain 
this evolutionary scenario by the following binary interac- 
tion. When the primary evolves and fills its Roche lobe, 
accreting matter to the secondary, if the accretion time- 
scale, = M/M, is shorter than the thermal time-scale 
of the secondary, tkh = GAP /{RL), both components can 
fill their Roche lobes, resulting in a contact binary and 
L2 Roche-lobe overfiow. The material passing through the 
L2 Lagrangian point, will remove angular momentum from 
the system, forming a disk-like structure around what will 
become a merger remnant of an interacting binary. This 
would imply that VFTS698 and R4 are in respectively a 
pre- and post-merger evolutionary state. However the al- 
ternative scenario sees the loss of angular momentum slow 
the system down enough to prevent a binary merger as the 
initially more massive companion will explode as a type 
Ib/c supernova. The result sees the mass donor slim down 
to a Wolf-R ayet star, lea. v ing a short-period binary system 
as found by iGlark et al.l (1201 1[) in Weste rlund 1 (further 
details can be found in iRitchie etaI1[20Tl . 

In Sect. 14.31 the nitrogen abundances were investigated 
assuming each component to be that of a stellar photo- 
sphere. However, in view of the above scenario, it is dif- 
ficult to ascertain what constraints can be placed on the 
analysis of nitrogen abundances. For the hot component. 



^ a possible ex c eption being the SMC star N 82, which 
iHevdari-Malaveril (| 19901 1 determined to be a sgB[e], but 
this classification was w eakend by the investigation of 
IMennickent fc SmithI (|2010l ) from an interacting binary nature. 
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our nitrogen abundance estimate of 8.0 dex still appears 
to be a reasonable result given the assumptions made for 
surface gravity and microturbulence. Considering the cool 
component is thought to originate from an optically thick 
disc, our assumption of a supergiant surface gravity can no 
longer be applied, and hence it is unclear what conclusions 
we can draw on the parameters of the pseudo-photosphere 
other than there is evidence for a temperature consistent 
with a mid-B type photosphere, with spectral evidence for 
a large nitrogen enrichment. 

6. Conclusions 

A detailed spectroscopic and photometric study of 
VFTS698 has been presented to ascertain its evolution- 
ary status. In summary, VFTS698, shows evidence for be- 
ing a double-lined spectroscopic binary, consisting of early- 
and mid-B type components. The spectrum is rich with 
iron-group emission and absorption features which appear 
syno nymous with the B[e] phenomenon (e.g.. lLamers et aLl 
11998) . OGLE-III photometry shows a stable 12.7 d period 
from 8 years of observations, and IRSF and Spiteer imaging 
shows VFTS698 to have strong near- and mid-IR excesses. 
We present the scenario for a W Serpentis-type interacting 
binary based on the following attributes: 

1 The near- and mid-IR excess of VFTS698 appears dif- 
ferent to that observed in bona fide sgB[e] stars; 

2 The stable 12.7 d photometric period is consistent with 
both the orbital motion and short-term orbital periods 
of W Serpentis-like objects. RVs of both components 
appear consistent with this short-period orbit; 

3 There is evidence for long-term photometric varia- 
tions, a definin g characteristic for this class of binaries 
(|Strupatlll987[ ): 

4 A complex circumstellar structure is inferred from the 
time- variable emission lines which display diverse mor- 
phologies, e.g., P Cygni Hei profiles, and shell-like Fell 
lines, similar to those in /? Lyrae (see|Tarasov 2000). 

The following points require contemporaneous photom- 
etry and spectroscopy, together with detailed modelling in 
order to better ascertain the exact nature of the VFTS698 
system. 

1 The near- and mid-IR excess, iron-group shell spectra, 
and P Cygni structure in the helium spectra are indi- 
cators for a complex circumstellar geometry consisting 
of circumbinary material, a circumstellar disc and wind 
outflows. 

2 The RV measurements of the primary /circumstellar 
disc appear inconsistent with an anti-correlation which 
would be expected for a binary system. 

Through the course of this investigation, a number of 
different comparison objects (see Table [7]) have been used 
to conclude that VFTS698 is an interacting binary, com- 
prising a visible, hot secondary in orbit with a veiled, cool, 
massive primary. We surmise the spectral component asso- 
ciated with the primary to originate from the circumstellar 
disc, but require more detailed RV measurements and con- 
temporary photometry to support this argument. 
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Table 7. Comparison of orbital and stellar parameters of VFTS698 derived here with those for known sgB[e] stars and 
interacting binary systems. 
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Table 8. Line identifications for the three plasma components of the 
VFTS698 system. Equivalent width (EW) and full-width half-maxinnim 
(FWHM) measurements were found by Gaussian-profile fitting, and 
are taken from the first epoch of observations with the LR02 setting. 
Identifications marked by an astrerix have insufficient evidence for an 
accurate identification. 
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0.6 


0.1 


4836.2 


Cm 


63.2 


5.0 


0.7 


0.1 


4848.2 


Cm 


101.5 


8.7 


0.6 


0.1 


4871.6 


Fell 


55.1 


11.0 


0.9 


0.1 


4877.5 


Cm 


140.0 


9.5 


0.7 


0.1 


4889.6 


[Fell] 










4894.9 


Cm 


18.2 


6.7 


0.5 


0.1 


4901.7 


Cm 


46.0 


9.6 


1.0 


0.2 


4905.4 


[Fell] 










4926.7 


Cm 










4952.7 


Fell 


28.7 


3.9 


0.7 


0.1 


4954.0 


Fell 


33.4 


3.6 


0.7 


0.1 


4958.2 


[Fen] 










5006.7 


[Fen] 










5019.5 


Fen 


240.9 


9.5 


0.6 


0.1 


5020.2 


[Fen] 










6446.4 


Fen 


107.2 


6.4 


0.6 


0.1 


6456.4 


Fen 


306.8 


6.2 


0.6 


0.1 


6482.2 


Fen 


106.3 


5.8 


0.5 


0.1 


6487.3 


Fen 


72.3 


6.9 


0.7 


0.1 


6614.3 


Fen 


127.3 


2.4 


0.6 


0.1 


6716.2 


Tin 


140.9 


2.5 


0.6 


0.1 


6731.3 


Fen 


95.0 


2.4 


0.6 


0.1 


6752.0 


Cm 


12.5 


2.1 


0.5 


0.1 


6764,1 


Fell 


22,9 


2,8 


0.8 


0,1 


()799.9 


i-'c 11 


12.5 


2.3 


0.() 


0.1 


Lines associated with the cool component: 


3995.00 


Nn 


198.7 


15.5 


1.3 


0.1 


4009.26 


He I 


425.6 


46.1 


1.1 


0.1 


4026.19 


He I 


523.7 


77.7 


1.1 


0.1 


4106 


Nil* 


166.5 


16.7 


1.3 


0.1 


4120.82 


He I 


444.2 


43.9 


1.0 


0.1 


4128.05 


Sin 


123.9 


12.8 


1.2 


0.1 


4130.89 


Sin 


122.7 


10.8 


0.8 


0.1 


4143.76 


He I 


421.8 


50.4 


1.0 


0.0 


4168.967 


He I 


187.9 


12.6 


1.0 


0.1 


4387.93 


He I 


555.6 


24.7 


1.1 


0.0 


4395.94 


On* 


217.9 


15.1 


1.3 


0.1 


4419.6 


Fein 


188.5 


9.9 


0.9 


0.0 


4431.02 


Fein 


57.6 


11.4 


0.7 


0.1 


4437.551 


He I 


287.3 


4.3 


1.1 


0.1 


4471.47 


He I 










4481.23 


Mg n 


93.7 


7.5 


1.1 


0.1 


4552 


Sim 


116.9 


14.5 


1.4 


0.1 


4567.84 


Sim 


157.2 


48.7 


2.2 


0.6 


4601.48 


Nn 


107.7 


10.6 


1.8 


0.1 


4607.15 


Nn 


118.2 


15.1 


1.9 


0.2 


4613.87 


Nn 


54.8 


10.0 


1.6 


0.2 


4643.09 


Nn 


275.7 


15.9 


2.6 


0.1 



continued on next page 



p. R. Dunstall et al.: A peculiar B[e]-like supergiant in 30 Doradus 
Table 8. continued 



Wavelength Species EW A(EW) FWHM A(FWHM) 



(A) 




(mA) 


(mA) 


(A) 


(A) 


4713.16 


He I 










4921.93 


He I 










Lines associated with the hot component: 


4008.00 


He I 


282.2 


44.8 


1.4 


0.2 


4088.84 


Si IV 


230.8 


12.6 


1.6 


0.1 


4097.33 


Nm 


141.5 


8.7 


1.3 


0.1 


4116.1 


Si IV 


135.1 


12.2 


1.4 


0.1 


4119 


He 


116.5 


20.3 


0.8 


0.1 


4143 


He I 


368.8 


48.2 


1.5 


0.2 


4167 


He I 


58.8 


11.6 


0.9 


0.2 


4387 


He I 


164.7 


21.3 


0.9 


0.1 


4437 


He I 


91.3 


20.7 


1.1 


0.2 


4567.84 


Sim 


104.0 


25.6 


1.3 


0.2 


4601.48 


Nil 


28.8 


9.9 


1.1 


0.3 


4607.15 


Nil 


24.0 


7.3 


0.9 


0.2 


4643.09 


Nil 


23.6 


9.3 


1.0 


0.3 


4685.71 


Hen 


47.3 


7.9 


2.0 


0.3 


4713.16 


He I 










4921.93 


He I 











